Edge spin excitations and reconstructions of integer quantum Hall liquids 
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We study the effect of electron-electron interaction on the charge and spin structures at the 
edge of integer quantum Hall liquids, under three different kinds of confining potentials. Our 
exact diagonalization calculation for small systems indicates that the low energy excitations of 
v = 1 ferromagnetic state are bosonic edge spin waves. Instabilities of the ferromagnetic state with 
altering confinement strength result from the softening of these edge spin waves, and formation of 
edge spin textures. In v < 2 regime, exact diagonalization on edge electron systems indicates that 
compact Hartree-Fock states with different total spin always become ground states in some regions 
of parameter space, and the ground states appear in between two compact states are their edge spin 
waves. The initial v — 2 instability is toward the compact state with total spin 1. Larger systems 
are studied using a microscopic trial wave functions, and some quantitative predictions on the edge 
instabilities for a certain type of confining potential are reached in the thermodynamic limit. 



I. INTRODUCTION 

In quantum Hall (QH) systems the charged bulk exci- 
tations are gapped, while gapless excitations exist at the 
edgesJ^ Due to the quantization of kinetic energy un- 
der strong magnetic field, the physics of edge excitations 
depends on the interplay of electron-electron interaction 
and the confining potential. The different manners in 
which the two dimensional electron gas (2DEG) is con- 
fined lead to rich electronic structures at the edges of the 
sample. 

The edge of a v = 1 spinless QH liquid confined 
in a simple geometry by a sharp confining potential is 
described by the Zp = 1 chiral Fermi-liquid theory^ 
With smoother confinement, the edge undergoes a charge 
reconstruction transition in which separated electronic 
lumps appear and bring new boundaries with even num- 
ber, breaking the chiralityj^ Further including the spin 
freedom and with sufficiently weak Zeeman splitting, 
Hartree-Fock (HF) technique shows that the charge re- 
construction transition is pre-emptied by formation of 
edge spin textured In fact closely related edge spin re- 
construction was found for v = 2 within HF approxi- 
mation, in which a spin polarized strip is created along 
the edge and through which the edge undergoes a second- 
order spin-unpolarized to spin-polarized transition as the 
confining potential smoothes.-^ The edge modes associ- 
ated with this spin reconstruction of v = 2 edge show 
up in NMR 8 and momentum-resolved electron tunneling 
spectroscopy^ Still within HF approximation, Barlas et 
al£ found that under sufficiently weak triangular confin- 
ing potential, unpolarized v = 2 state's spin-polarization 
transition is preceded by a charge reconstruction of a sin- 
gle spin species, creating a detached spin-polarized strip. 
This kind of edge structure are understood as a simul- 
taneous edge charge and spin reconstruction. It is these 
recent theoretical and in particular, experimental devel- 
opments that motivated the present work, which go be- 
yond earlier theoretical work^i^r— in several aspects. 

First of all, earlier studies of integer QH edge instabil- 



ities and reconstructions were mostly based on HF ap- 
proximation, and focused on the ground states only. In 
the present work we perform detailed exact diagonaliza- 
tion study on the same systems. In addition to obtaining 
the exact ground states, we also obtain the low-energy 
spectra of the system before instabilities occur, which 
provide deeper insight into the nature of the instabilities 
and the mechanisms of the corresponding edge recon- 
structions. 

Secondly, details of a QH liquid's edge structure may 
also depend on the detailed form of the confining poten- 
tial. Earlier works usually focus on one specific type of 
confining potential. In this paper, we consider a 2DEG 
with disk geometry under three different types of con- 
finement. Positive background charge confining poten- 
tial is usually used to model a realistic confinement. 13 
Parabolic confining potential is commonly used in the 
research of QH dots4^ By numerically solving the Pois- 
son and Schrodinger equations for electron states in 
a GaAs/Al x Ga\- x As heterostructure with three spa- 
tial dimensions' confinement in Hartree approximation, 
Kumar et al. argue that it is reasonable to model a 
real QH dot device's confining potential as a parabolic 
oneJ£ Linear confining potential also deserves some at- 
tention, because the recent momentum-resolved tunnel- 
ing experiments for parallel wires which can be modeled 
as linear confining potential, have observed edge spin 
reconstruction, and more complicated edge charge and 
spin reconstruction is found in a similar triangular con- 
fining potential^ Such a comprehensive study allows us 
to compare the systems' behavior under different types 
of confinements, and distinguish between universal and 
non-universal aspects of edge physics. 

Last but not least, here we study v = 1 and v = 2 
edges together, and reveal the similarities and differences 
in their behavior. Such side-by-side comparison has not 
been performed before. In fact insights from the study of 
v = 1 are crucial to our understanding of the v = 2 case. 

Our most robust results are summarized as what fol- 
lows. Ignoring Zeeman coupling, the low energy spectra 
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of v = 1 ferromagnetic state can be mapped onto those 
of AS* = — 1 bosons on top of the ferromagnetic state. 
These excited bosons are edge spin waves (ESWs) , which 
propagate in both directions. Under positive background 
charge confining potential, the ferromagnetic state is sta- 
ble as the distance d from positive background charge 
layer to 2DEG approaches zero. When the distance d be- 
comes larger, ferromagnetic state is destabilized by soft- 
ening of ESWs, resulting in edge spin textures. Under 
parabolic or linear confinement, the ferromagnetic state 
is destabilized by softening of ESWs under both smoother 
and stronger confinements, and the ferromagnetic state 
window disappears as the particle number iV increases. 
In v < 2 regime, compact states (defined below) with 
different total spin S can always become ground states 
in some regions of parameter space. The low energy ex- 
citations of a compact state can be mapped onto those of 
AS = — 1 bosons on top of it, thus are ESWs of this com- 
pact state. The ground states appearing in between two 
different compact states can thus be viewed as edge spin 
textures as well. The initial v = 2 instability is toward 
the S = 1 compact state. 

The rest of the paper is organized as what follows. 
In Sec. II, we introduce the models, provide some nu- 
merical details of our calculations, and an electrostatic 
model which gives some qualitative understanding on in- 
teger quantum Hall liquid's edge instabilities. In Sec. 
Ill we investigate the low energy excitations and edge 
instabilities of v — 1 ferromagnetic state using exact di- 
agonalization on small systems and a microscopic trial 
wave function on large systems. Sec. IV considers the 
low energy excitations and edge instabilities of v = 2 
unpolarized state. Besides carrying out exact diagonal- 
ization on small systems, we also study the large systems 
by separating out the edge electron systems and using ex- 
act diagonalization on them. Some concluding remarks 
are made in Sec. V. We make comparisons with earlier 
work wherever appropriate. 



II. THE MODELS AND ELECTROSTATIC 
CONSIDERATION 

A. The models 

We consider a 2DEG with spin degree of freedom in a 
disk geometry. The single particle Hamiltonian is 
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-A(r)} 2 + U(v)+ gf i B s z B, (1) 



where m e , — e is the mass and charge of a single elec- 
tron respectively, c is the speed of light, /is is the Bohr 
magneton, p is the momentum operator. U(r) is the ro- 
tationally invariant confining potential. The vector po- 
tential A = (—By/2,Bx/2,0) in symmetric gauge. s z 
is the z axis component of single particle spin operator. 
Confining electrons to the lowest Landau level (LL), the 



single particle wave functions are 

cj> m (z) = (27Tl 2 B 2 m m\r 1 / 2 z m e-^ 2 ^, (2) 

where z — (x — iy)/ls is the complex coordinate in the 
plane of the 2DEG and Ib — \J hc/eB is magnetic length, 
to = 0, 1, 2, • • • . The complete Hamiltonian in symmetric 
gauge is then 
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where cj^ CT is the electron creation operator for the lowest 
LL single-electron state with orbital angular momentum 



to and spin er, h r , 



is the occupation number 



operator of the TOth orbital with spin a, and 
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J J eri2 

(4) 

represents the electron-electron Coulomb interaction (e is 
the dielectric constant). C/^f is the matrix element of the 
rotationally invariant confining potential. Hamiltonian 
(3) has rotational symmetry, and also commutes with 
total spin operator S and its z axis component S z . Total 
angular momentum M, total spin S and S z are good 
quantum numbers. 

In GaAs, arising from band structure and spin-orbit 
coupling, electron spin g factor is renormalized to g 
—0.44. The ratio of Zeeman energy gusB to the typi- 
cal Coulomb interaction energy e 2 /(eZs) is nearly 1/55 
under magnetic field B = 9T. Thus the Coulomb inter- 
action dominates the Zeeman gap in the magnetic fields 
of interest. Actually, the g factor could even be tuned by 
quantum confinement^ or by application of hydrostatic 
pressure, and could pass through zero in these circum- 
stances. We thus ignore the Zeeman term when studying 
the edge spin excitations of integer quantum Hall liquids 
most of the time. As a result for each energy level with 
quantum number S obtained in exact diagonalization, 
it has degeneracy 2S + 1 with different choices of S z . 
Adding back a Zeeman term will not change the eigen- 
states, but only split the original degenerate energy lev- 
els and give the states with different S z different energy 
shifts; we will also consider the effects of such splitting 
below. 

To model the confinement of 2DEG based on a 
modulation-doped AlGaAs/GaAs heterostructure, we as- 
sume a uniformly distributed neutralizing positive back- 
ground charge layer at a distance d above the 2DEG 
as in Fig. 1. The radius of the positive background 
is R = y/2N/v, so the disk encloses exactly N/v mag- 
netic flux quanta for iV-particle 2DEG system. Coulomb 
interaction between the positive background charge and 
the 2DEG gives rise to the positive background charge 
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FIG. 1: (Color online) Sketch of an electron (blue) layer and 
a uniformly distributed neutralizing background charge (red) 
layer separated by a distance d from each other in disk geom- 
etry. 
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where p is the charge density of the background. The 
ratio of distance d between the electron layer and the 
background charge layer to magnetic length lg: d/ls, is 
the tuning parameter that controls the relative strength 
of the confining potential to the electron-electron inter- 
action. 

For parabolic confining potential U (r) = f3r 2 , 



(m\ /3r 2 \m) n r , 



= bJ2 2(m + l)(^-) 
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The coefficient (3 represents the steepness of the parabolic 
confinement, which is in units of meV/nm 2 . In this case, 



h = pr B /(- r ) 

elB 
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is the dimensionless tuning parameter, which is propor- 
tional to steepness coefficient (3 and the 3/2 power of 
magnetic field strength B. The one body confining po- 
tential operator can be written as the function of total 
orbital angular momentum operator M, and total parti- 
cle number operator N . 
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where the coefficient a represents the slope of the linear 
confining potential, which is in units of meV/nm. Here, 
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is the dimensionless tuning parameter, which is propor- 
tional to the confining potential's slope a and magnetic 
field strength B. 

For the exact diagonalization calculations of small sys- 
tems below, we limit N particles in a finite number N or t, 
of orbitals, and solve the problem in subspaces with cer- 
tain quantum numbers M, S and S z . Therefore the 
Hilbert space could be limited to a reasonable size and 
exact diagonalization could be carried out. For the ex- 
act diagonalization calculations in Sec. III. A and IV. A, 
we give each electron system a sufficiently large orbital 
number N orb (N orb = 16, 23, 19 for N = 8, 10, 12 system 
respectively), so that the calculated orbital occupation 
numbers n m approach in the outermost given orbitals. 
For this reason, we believe that the limited orbital num- 
ber N or b we have chosen has little effect to restrict the 
electrons to the inner orbitals, and the results will not 
change qualitatively if larger N orb is used in this exact 
diagonalization calculation. 

At this point, we also need to introduce a standard to 
identify the filling factor v of the electron systems for the 
future calculation. For finite systems it is very easy to 
identify the v — 1 and 2 states, because electrons just 
completely occupy the innermost lowest LL orbitals of 
the single or two spin component (s) and form a Slater 
determinant state. The finite size systems for v = 2 
state study should have even numbers of particle. As 
the confinement smoothes, electron system transits from 
unpolarizcd v = 2 state to partially polarized states, and 
finally to polarized v — 1 state. When the system is in 
this transition region, we say it is in 1 < v < 2 regime. 
For an electron system with filling factor v = 1, when 
the confinement becomes smoother (stronger), electrons 
would tend to distribute in more (less) orbitals. If the 
confinement becomes a little smoother (stronger) and the 
system has just undergone a few transitions from the 
original v — 1 state, we say this system is in v < 1 {y > 1) 
regime. The identification is similar for v = 2 state. But 
filling factor could not be larger than 2 if Hilbert space is 
truncated to the lowest LL, so only v < 2 regime exists. 



B. Electrostatic consideration of origin of edge 
instabilities 



Before studying edge instabilities of v = 1 and v = 2 
QH liquids numerically, we introduce an electrostatic 
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FIG. 2: (Color online) (a) Edge test electron and electron 
layer in disk geometry, (b) Edge test electron and eletron layer 
in strip geometry, (c) Edge test electron and semi-infinite 
electron layer and positive background charge layer. 



model, which gives us some qualitative understandings 
of the instabilities, especially for large systems. The ex- 
pectation value of r 2 for single electron state of the low- 
est LL is 2(m + 1)1%, so it is reasonable to make the 
approximation that every 1 (2) electron(s) occupy area 
2irl B and the whole —Ne charge uniformly distributes in 
a disk with radius R = V2Nl B {VNI b ) for v = 1 (2) 
state. The charge density is a — — e/ '(2ttI b ) for v = 1 
state, and —e/(irl B ) for v = 2 state. 

The edge instability and the resultant structure is de- 
termined by the competition between electron-electron 
repulsion and confinement potential. We thus consider 
the forces felt by a test electron at the border of this 
disk. Quantum mechanically the electron's wave function 
has width ~ Is in the radial direction. For this reason, 
in the electrostatic model, we put the test electron at a 
point which is Iqm away from the electron layer's bor- 
der and Iqm ~ Ib- Besides direct interaction, exchange 
interaction between two electrons with same spin in the 
2DEG, is also an essential factor of edge instability. To 
include exchange interaction in our electrostatic model, 
we consider the following. First, exchange interaction is 
attractive, which can be modeled as the direct interac- 
tion between some hypothetical positive charges in the 
disk and the edge test electron. Second, exchange in- 
teraction has a short range with the order of magnetic 
length Ib in the QH regime. Therefore, including the 
effect of exchange interaction at v = 1 (where exchange 
effects exists between all pairs of electrons), we reduce 
the radius of the (negatively charged) disk from R to 



R— l ex without changing the charge density cr, where l ex 
also has the order of Ib- A similar model (with proper 
modification) can be used to describe the v — 2 case. 

As shown in Fig. 2(a), the electrostatic model for 
2DEG v = 1 state is an electron disk with uniform charge 
density a — —e/(2wl%), radius R — \/2NIb — l e x, an d an 
edge test electron with distance / = Iqm + lex away from 
the electron disk. For v = 2 state, we need two electron 
disks and one edge test electron to model the system; one 
of the disks with uniform charge density a — —e/(2nl%) 
and radius R — >/NIb ~ I ex, models the electrons in the 
plane with the same spin as the edge electron, and the 
other disk with uniform charge density a = —e/{2-Kl%) 



and radius R = \/NIb, models the electrons in the plane 
with the different spin as the edge electron. Edge test 
electron is I = l ex + Iqm away from the first disk and 
Iqm away from the second disk in the v = 2 state's elec- 
trostatic model. 




FIG. 3: (Color online) Numerically calculated force F felt by 
the test electron vs R/l in the range R/l £ [0, 500] (blue line) 
and its fitting function 2.0251 \n(R/l) (dashed red line) in disk 
geometry. 

To study the edge instabilities, we first consider the 
force felt by the edge test electron from the electron disks 
in this model. In the strip geometry of the same prob- 
lem, it is easy to obtain an explicit form of this force, 
which is what we consider first. As Fig. 2(b), it is 
straightforward to show that the edge test electron feels a 
force F = 2 \n.(L/V)(ea/e) in x direction from the electron 
layer, where L is the width of the strip and / is the dis- 
tance from the test electron to the electron layer. When 
the system is large (R^> Ib), the force felt by edge test 
electron in disk geometry tends to be the same as the 
one in strip geometry, where R plays the role of L. In 
disk geometry as Fig. 2(a), force F felt by the edge test 
electron vs R/l in range [0, 500] is calculated numerically 
in Fig. 3. It is well fitted by a function 2.0251 ln(R/l). 
The coefficient of \nR/l approaches 2 as the strip one 
with increasing R. Therefore, numerically we can see the 
force felt by edge test electron in disk geometry has a 
logarithmic dependence of R/l. I equals Iqm + lex as the 
test electron's spin is the same as the ones in electron 
disk, and Iqm as they are different. 

The inward force from the confining potential is pro- 
portional to R under parabolic confining potential, and 
is a constant under linear confining potential. Because 
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of logarithmical R dependence of the outward force from 
the electron disk, thermodynamic limit does not exist for 
edge instability critical points (b c and a c ) for these two 
types of confinements. By further approximating that 
Iqm ~ Ib, lex ~ and F ps 21n(iZ//), the electrostatic 
model predicts that for N ^> 1 the critical point of v = 1 
instability (smoother confinement side) decreases (in- 
creases) as b c ~ [ln(2iV)]/(47rV2A0 (o c - [ln(2JV)]/(27r)) 
with increasing particle number iV under parabolic (lin- 
ear) confining potential. By the same approximations for 
N 3> 0, the critical point of i/ = 2 instability decreases 
(increases) as b c ~ []nN]/(2ny/N) (a c ~ (hxN)/w) with 
increasing particle number N under parabolic (linear) 
confining potential. 

For large systems under positive background charge 
confining potential, the electron layer and background 
charge layer could be approximated as two semi-infinite 
layers as Fig. 2(c). For v = 1 state, when the distance of 
the two layers d — \Jl\ x + 21 £X Iqm, the force felt by the 
edge test electron becomes zero and v = 1 edge instabil- 
ity happens when d further increases. For v = 2 state, 
when the distance of the two layers d — y/l^ x /^ + lexlQM, 
the force felt by the edge test electron becomes zero and 
v = 2 edge instability happens when d further increases. 
Because Iqm and l ex are independent of R when i? 3> Ib, 
the thermodynamic limit exists for positive background 
charge confinement. From the electrostatic model, we 
can also predict that in thermodynamic limit, the v = 2 
edge instability critical point d c is larger than the one for 
v = 1 edge instability, and they are both of order I b ■ 

These conclusions are confirmed by numerical results 
presented below, and also consistent with earlier work on 
fractional QH edge reconstructions^ 
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FIG. 4: (Color online) Low energy spectrum of 10-electron 
system under parabolic confining potential when 6 = 0.06. 
The different total spin quantum numbers S of the eigenstates 
are labeled by different colors as the annotation. The ground 
state is v = 1 ferromagnetic state (F). Low energy excitations 
in AM > subspaces are enclosed by dashed boxes, and 
low energy excitations in AM < subspaces are enclosed by 
dotted boxes. The number below each box is the number of 
states inside the box. 
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III. LOW ENERGY EDGE EXCITATIONS AND 
EDGE INSTABILITIES OF v = 1 
FERROMAGNETIC STATE 

In this section, we investigate the low energy excita- 
tions of v — f ferromagnetic state under the three kinds 
of confinements. First exact diagonalization results are 
presented for small systems. The low energy spectra can 
be mapped onto those of AS = — 1 bosons on top of the 
ferromagnetic states. By investigating the spin structure 
of these low energy excitations, we conclude that they are 
edge spin waves (ESWs). Microscopic trial wave func- 
tions are constructed to describe the ESWs, which allow 
for their studies in larger systems. We conclude that the 
edge instabilities of the v = 1 ferromagnetic state are 
triggered by the condensation of the ESWs which results 
in edge spin textures. 



A. Exact diagonalization study for small systems 

In this subsection exact diagonalization method is used 
to study the small electron systems close to v = 1 fer- 



TABLE I: One-to-one correspondence between the configura- 
tions of spin AS = — 1 bosons and the low energy excitations 
in each angular momentum subspace AM of TV = 10 system 
under parabolic confinement. AM is the difference of angu- 
lar momentum quantum number M compared to the one for 



ferromagnetic state, 
momentum q. 



is the number of bosons with angular 



romagnetic state under the three confinements. Fig. 4 
is the fO-clcctron system's low energy spectrum under 
parabolic confinement with b = 0.21 (ferromagnetic state 
is the ground state). We will use this specific spectrum 
to show some low energy spectra's common character- 
istics shared by all small systems we studied and un- 
der all the three confinements. AM = M — Mq is the 
difference of orbital angular momentum quantum num- 
ber M compared to the one M = Mq for ferromag- 
netic state (Mo = 45 for N = 10 system). In the sub- 
spaces with AM > 0, a set of low energy excitations 
(enclosed by dashed boxes in Fig. 4) are separated from 
higher energy states. One of these excitations becomes 
the ground state and thus destabilizes the ferromagnetic 
ground state with smoother confinement. These low en- 
ergy excitations all have nonzero AS and the polarized 
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FIG. 5: (Color online) Charge density n c (red lines) and 
change of spin density As z (dashed red lines) for the low 
energy excitations of ferromagnetic state, in (a) subspace 
AM = 1, S z = S = 4 and (b) subspace AM = -1, 
S z = S — 4, of the 10-electron system under parabolic con- 
finement with b = 0.06. The charge density profiles of low 
energy excitations have no big change compared to the one of 
ferromagnetic state (blue line). 



excitations with AS = have higher energies compared 
to them. In the subspaces with AM < 0, another set of 
low energy excitations (enclosed by dotted boxes in Fig. 
4) are separated from higher energy states, and one of 
these excitations becomes the ground state with stronger 
confinement (except for positive background charge con- 
finement case where the ferromagnetic state remains sta- 
ble all the way to d — 0.). Polarized excitation with 
AS = does not exist in the subspaces with AM < 0, 
because all the lowest orbitals with up spin are already 
occupied in the ferromagnetic state. More interestingly, 
the low energy spectrum (enclosed by the boxes) matches 
that of a bosonic system, with each boson reduces the 
system's total spin by 1. The one-to-one correspondence 
between the bosonic states and low energy excitations of 
this 10-electron system under parabolic confinement is 
shown in Table I. Figs. 5 shows the charge and spin den- 
sity profiles of the bosonic states in subspaces AM = 1 
and —1 of the spectrum Fig. 4. While charge density 
profiles are almost identical to that of the ferromagnetic 
state, we find the flipped spin is confined to the edge. We 
thus conclude the bosons excited from the ferromagnetic 
state are actually edge spin waves (ESWs)ri which can 
propagate both along the clockwise (AM > bosons) 
and anti-clockwise (AM < bosons) directions. We re- 
peat that the characteristics stated above are shared by 
all small systems we studied and under all the three con- 
finements. 

Fig. 6 presents the phase diagrams of 8, 10-electron 
systems close to v = 1 ferromagnetic states under 
the three confinements. As the confinement becomes 
smoother, the energy of a single boson state with some 
positive angular momentum q — qo corresponding to 
clockwise ESW first crosses that of the ferromagnetic 
state, thus destabilizes it. There are interactions among 
these bosons. As the confinement further smoothes, the 
bosons start to condense because in this process the 
benefit of adding a new boson with angular momentum 
qo would exceed the cost of the accompanying increase 
of interaction energy. When the confinement becomes 
smoother, the minimum of boson dispersion shift from 



q = qo to q = q' . If that happens, the system will 
choose to condense q' bosons, instead of go bosons. Us- 
ing this picture, we can understand why in the phase di- 
agrams Fig. 6(a-c) the systems flip spins one by one with 
smoother confinement and quantum number AM + qoAS 
is conserved for all the ground states, qo may change 
along with different tuning parameters as in Fig. 6(c). 
When the confinement becomes stronger, the systems 
start to condense bosons with q = — 1, which correspond 
to anti-clockwise ESWs. As in Figs. 6(d,e), the systems 
will condense all the possible q = — 1 bosons to form spin 
singlet states either at one time, or within a much smaller 
range of tuning parameter compared to the transitions in 
smoother confinement side. 

Exact diagonalization results also show that under pos- 
itive background charge confining potential, once the sys- 
tem transits to v = 1 ferromagnetic state, the phase will 
be stable as d/ls approaches zero. 



B. Microscopic trial wave function study for large 
systems 

Having identified the low energy excitations at v = 1 
as bosonic ESWs on top of the ferromagnetic state, we 
want to construct microscopic trial wave functions to 
approximate them. Like Abelian bosonization in one- 
dimensional Fermi gas, the bosonic generators do not re- 
spect the SU(2) spin rotational symmetry. For this rea- 
son, we only choose one state among all the degenerate 
states in each energy level with spin quantum number 
S to represent all states in this level. The state we will 
choose is the one with S z = S, because upon adding a fi- 
nite Zeeman term this state's energy is the lowest within 
this family. 

On top of the S z = S ferromagnetic state 
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the bosonic generators of ESWs were first constructed by 
Oaknin et alJ£ We rewrite them as 



N-l 



(11) 



m=0 



where b} n q cr = c f m+q (7 c m ^, q is an integer and a =f or 
4,. The coefficients ifj m , q ,a are obtained by minimizing the 
energy of B^ a \F), and the condition that it is normal- 
ized. There is one spin flip if a = j. (AS Z = — 1), and no 
spin flip if a =f (AS Z — 0). Because 



[Bq >a , 13^, a ,] — 3q,q' &a,(r' ^ ^ ^m,q,a (^m.q.q- ^m,g.q ^m,q,er ^m,g,(7 ) 
m=0 

(12) 



but 



(F\[B q , a ,Bl, !<T ,}\F)=6 g , g ,6 <T , a , 



(13) 
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FIG. 6: (Color online) Phase diagrams of 8, 10-electron systems under (a) v = 1 positive background charge, (b) parabolic, (c) 
linear confining potential in v < I regime; and under (d) parabolic, (e) linear confining potential in v > f regime based on exact 
diagonalization. AM is the difference between orbital angular momentum quantum number M and that of v = I ferromagnetic 
state for each system. The different total spin quantum numbers S of the eigenstates are labeled by different colors as the 
annotation. Step values are A(d/?s) = I0~ 2 , Ab = 2 x f0 -4 and Aa = 2 x f0~ 3 for the three confinements respectively. The 
configuration of bosons is written beside each phase and n q is the occupation number of the boson with angular momentum q. 



the states generated by i3| a can be approximated as 
bosons, especially when the number of excitations are 
small. The state of n bosons with q and a is 



N-l 



(St ff )« |F) = W^cL +g , CT c m , t )" \F) 



m=0 



(14) 



Compared with the exact low energy excitations, the 
microscopic trial wave functions (|14[) have good quantum 
numbers M and S z but not S. Nevertheless they do have 
large overlaps with the exact excited states (see Table 
II for examples). The overlaps increase with increasing 
particle number N, and we do not fully understand the 
reason of this surprising behavior yet. 

The bosonic excitations (|T4j) are classified into two dif- 
ferent types. The excitations generated by ^ (AS Z = 
sector) are the well studied edge magneto-plasmons 
(EMPs)£4 Since B^ t \F) = for q < 0, EMP can only 
carry positive angular momentum and is thus chiral. It 
excites the edge by creating a charge density wave with 
positive angular momentum q which propagates along 
the clockwise direction only. The excitations generated 
by ^ (AS Z — —1 sector) are the ESWs discussed in 
last subsection. As a non-chiral excitation, ESW excites 
the edge by creating a clockwise edge spin wave with 
positive q or anti-clockwise edge spin wave with negative 

q- 
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positive background 


parabolic 


linear 
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0.8643 


0.8584 


0.8865 
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0.9164 


0.9136 


0.9134 


10 


0.9889 


0.9896 


0.9530 


11 


0.9904 


0.9910 


0.9864 


12 


0.9916 


0.9920 


0.9875 


13 


0.9926 


0.9929 


0.9884 


14 


0.9934 


0.9936 


0.9891 


15 


0.9941 


0.9942 


0.9897 


16 


0.9946 


0.9947 


0.9901 


17 


0.9951 


0.9951 


0.9905 


18 


0.9955 


0.9954 


0.9909 


19 


0.9959 


0.9957 


0.9912 


20 


0.9962 


0.9960 


0.9914 



TABLE II: Overlaps between B\ ^ \F) and the exact lowest en- 
ergy state in the subspace with M — Mo + 1 and S z = N/2 — 1, 
as functions of the number of electrons under three different 
kinds of confinements. In calculating these states, tuning pa- 
rameters of confinement are chosen as the critical values of the 
ferromagnetic state's initial instibility with smoother confine- 
ment. 



For larger systems, exact diagonalization is not feasi- 
ble. The microscopic wave functions (|14|) give us a tool 
to study the instabilities near the ferromagnetic states. 
The energies of a single boson with q, a and two bosons 



8 



with the same q, a on top of the ferromagnetic state are 
(F\ B q , a HB\ ia \F) 



{F\B q , 



rBl* \F) 



-(F\H\F) (15) 



and 



AE, 



(?,°0 



{F\{B q ^H{Bj a f\F) 
(F\(B q , a f(Bl a y\F) 



(F\H\F) (16) 



respectively, e 



AE^'^ is the kinetic energy of a 



single boson, and e 



= AE^' a) -2AE] q '°' is the inter 



action energy between the two bosons. The energy of n 
bosons with the same q and a on top of the ferromagnetic 
state is 



ne 



(g,a) , n{n~ 1) ( g , g ) 



+ 



(17) 



In the exact diagonalization calculation for small sys- 
tems, besides ferromagnetic state, the ground states are 
always the bosonic states consisting of bosons with the 
same q and a. Assuming this is also true for larger sys- 
tems and using Eq. (17), we obtain the phase diagrams 
close to ferromagnetic states under the three confine- 
ments for particle number up to 40 (Fig. 7). We first no- 
tice that when the confinement smoothes, phase bound- 
aries of ferromagnetic states are determined by the ESW 
instabilities, rather than the EMP instabilities. With 
smoother confinement, ESWs are generated one by one 
to lower the systems' energy. 

For positive background charge confining potential, 
ESWs in AS Z = — 1 sector with q = 1 are generated 
on the ferromagnetic state one by one with smoother 
confinement. The critial value of ferromagnetic state's 
initial instability approaches a constant with increasing 
particle number, which implies the existence of thermo- 
dynamic limit. Therefore, we predict that in thermody- 
namic limit, the instabilities close to ferromagnetic state 
under positive background charge confinement are the 
condensation of AS Z = — 1 bosons with modulation an- 
gular momentum q = 1 (the lowest possible modulation 
angular momentum), and the critial point of ferromag- 
netic state's initial instability <1 c /Ib ~ 0.8. 

Under parabolic and linear confinements, the systems 
tend to condense AS Z = — 1 bosons with larger modula- 
tion angular momentum q to minimize the energy, espe- 
cially when the particle number becomes large. For a cer- 
tain system, the minimum of the bonson's dispersion may 
also change from one angular momentum q to another. 
Fig. 8 shows the dispersion relation of ESW in N = 40 
system under parabolic confinement when b = 0.0387. 
At this tuning parameter, the minimum is at q = 2 and 
its energy is lower than the ferromagnetic state's energy, 
so the system prefers to condense q — 2 bosons. With 
smaller b, the minimum would shift to q = 3, then the 
system gives up the q = 2 bosons and chooses to condense 
q = 3 bosons. This process can also be observed in the 
phase diagram Fig. 7(b), where the color of instabilities' 
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FIG. 7: (Color online) Phase diagrams based on trial wave 
functions close to v = 1 for large systems under (a) v = 1 
positive background charge, (b) parabolic and (c) linear con- 
fining potentials, v — 1 ferromagnetic states (F) transit to 
bosonic states of the form (B^ j_) n |_F) with varying strength 
of confinement. Whenever one more boson is generated, a 
point is plotted to label the critical tuning parameter of this 
transition. The color of each point represents the angular mo- 
mentum quantum number q of generated bosons, as labeled in 
the annotation. There is one exception. Under parabolic con- 
finement and with stronger confinement, the ferromagnetic 
state condenses q = — 1 bosons one by one when N < 12, 
but transits to unpolarized (even TV) or spin 1/2 (odd JV) 
state directly when N > 12 (its critical point is labeled by 
a bigger gray point). In the smoother confinement side, the 
phase boundaries of ferromagnetic states are determined by 
the edge spin wave instabilities, rather than the edge magne- 
toplasmon instabilities (labeled by triangles). For 8 and 10- 
electron systems, the critical points of ferromagnetic state's 
initial instability calculated by exact diagonalization are la- 
beled by crosses in smoother confinement side, and circles in 
stronger confinement side. 



critical points changes from red to green with smaller b 
when N — 40. As the confinement is strengthened un- 
der parabolic and linear confinement, AS Z — — 1 bosons 
with q = — 1 are always favored to be condensed. It is the 
condensation of these ESWs that drives spin edge recon- 
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AE(e z /d B ) 

0.12 




FIG. 8: Dispersion relation calculated by trial wave functions 
for edge magnetoplasmon (EMP) and edge spin wave (ESW) 
in N = 40 system under parabolic confining potential (b = 
0.0387). The arrows label the minimum of each branch of 
excitations. The dispersion relation of ESW is fitted well by 
a parabola. 



struction (or formation of edge spin textures) at v = l. 17 
First order phase transition will happen under parabolic 
confinement (N > 12). Instead of generating bosons one 
by one as the linear confinement case or small systems 
(N < 12), the system prefers to generate all the possible 
bosons at one time. The window of v = 1 ferromagnetic 
state decreases to zero with increasing particle numbers 
under parabolic and linear confinement. Thus the v = 1 
ferromagnetic state will not appear as ground state in 
large systems under these two confinements. Also con- 
sistent with the electrostatic analysis, the critical points 
of edge instabilities in the smoother confinements side 
decrease (increase) with increasing particle number un- 
der parabolic (linear) confining potential. For 8 and 10- 
electron systems, our trial wave function method's results 
are in quantitative agreement with exact diagonalization 
results of Sec. III. A; the latter are labeled in Fig. 7 for 
comparison. 
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FIG. 9: (a) Stability diagram calculated by trial wave func- 
tion for the 40-electron system's ferromagnetic state (F) un- 
der positive background charge confining potential, showing 
the region of stability to the ferromagnetic state, and the re- 
gions where it is unstable to edge spin wave (ESW) and edge 
magneto-plasmon (EMP) instabilities, (b) Renormalized crit- 
ical g factor § c = g c fiBB/(e 2 /eIb) for the ferromagnetic state 
having EMP, instead of ESW initial instability vs particle 
number N under positive background charge confining po- 
tential calculated by trial wave function. 



Finite Zeeman term will lift the energies of ESWs. The 
dimensionless parameter g — gugB / (e 2 / cIb) is the ratio 
of the Zeeman energy to the typical Coulomb energy. 
For large enough g — g c the ferromagnetic state will pre- 
fer to transit to its polarized excitations, and EMP will 
replace ESW to become the initial instability of ferro- 
magnetic state. Fig. 9(a) shows the stability diagram of 
the 40-electron system under positive background charge 
confinement in (g, d/ls) plane with phase boundaries 
obtained by our trial wave function. For small Zeeman 
energies, g < g c = 0.0056, the ferromagnetic state is 
destabilized by the ESWs with smoother confinement. 
As g c increases, the critical tuning parameter <1 c /Ib of 
ferromagnetic state's initial instability increases as well 
until at g c = 0.0056, the polarized EMP becomes the ini- 
tial instability. Fig. 9(b) shows how the critical value g c 
changes with increasing particle number. With thermo- 
dynamic limit under positive background charge confine- 
ment g c approaches a constant with increasing particle 
number. Based on this figure, we predict that in thermo- 
dynamic limit, g c m 0.005. And if we choose dielectric 
constant e = 12.8, g = -0.44 for the 2DEG based on 
a modulation-doped AlGaAs/GaAs heterostructure, the 
critical magnetic field magnitude is B c — 0.74T corre- 
sponding to the renormalized critical g factor g c = 0.005. 
Our quantitative prediction of g c is about one half that 
obtained by J. Sjostrand et al.^ who used a modified 
version of the positive background charge confinement 
which is slightly sharper. Other results like d c also differ 
by a order 1 numerical factors. These indicate such quan- 
titative results are sensitive to details of the confining po- 
tential. The fact that g c is larger for sharper confinement 
is not surprising as sharper confinement tends to lift the 
energy of charge modes more than spin modes. We note 
that in real samples there are usually additional sources 
of confining potential on top of the background charge, 
like gates and in particular, sharp crystalline boundaries. 
Thus in reality B c may well exceed ~ IT. Also in real- 
ity d <C d c . As a result edge spin textures (possibly in 
combination with charge reconstruction) are likely to be 
present in many systems. 



IV. LOW ENERGY EDGE EXCITATIONS AND 
EDGE INSTABILITIES OF v = 2 UNPOLARIZED 
STATE 

In this section, we investigate the low energy excita- 
tions and edge instabilities of v — 2 unpolarized state 
under the three kinds of confinements. Carrying out ex- 
act diagonalization on small systems in 1 < v < 2 regime, 
we observe that compact states (defined below) with total 
spin from to N/2 become ground states in some regions 
of tuning parameter. Similar to the ferromagnetic state 
case, low energy spectra can be mapped onto those of 
the bosons on top of the v = 2 unpolarized state. The 
non-compact ground states with simultaneous edge spin 
and charge reconstruction are understood as the conden- 
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sation of bosons (ESWs) on top of S ^ compact states. 
Larger systems are studied by separating edge electrons 
from bulk ones; this allows us to reach conclusions about 
thermodynamic limit in certain cases. 



A. Exact diagonalization study of small systems 

In this subsection, we study small systems in 1 < v < 2 
regime. Exact diagonalization is carried out in 8, 10, 12- 
electron systems, and the phase diagrams of 10-electron 
system (Fig. 10) are used to show some common char- 
acteristics shared by all small systems we studied. Start- 
ing from v = 2 unpolarized state (with M — 20 in Figs. 
10), with smoother confinement, the compact states with 
spin S from to N/2 become ground states (labeled by 
black points) in some tuning parameter regions. A com- 
pact state is the state having the minimum M compatible 
with a given S (the degeneracy is 2S + 1 without Zeeman 
coupling). The compact state with S z = S is 



N-k-1 



fc-i 



\c k N - k )= n 4 )t n 



|o> 



(18) 



m— 



m'=0 



It is a single Slater determinant with the N — k lowest 
m single particle states with spin up and the k lowest m 
single particle states with spin down occupied, in which 
k is an integer running from to N/2. They are the 
ground states obtained in previous work based on HF 
approximation^! The fact that such compact states oc- 
cupy large portions of the phase diagrams in our exact 
diagonalization study provides strong support to earlier 
HF studies. The compact states with S z < S can be 
obtained using spin lowering operator 



s- = 



W 

/ y m 



(19) 



on the S z = S compact state |C^_ fc ). Besides com- 
pact states, other ground states with smaller S (labeled 
by red points in Figs. 10) appear between neighbor- 
ing compact states. These non-compact ground states 
have charge and spin textures, and will be studied in the 
next subsection. We find that the difference between a 
non-compact ground state and the neighboring compact 
states \C N _ k ) and \C^~_} k+1 ) lies mostly in the occupa- 
tion numbers between the (fc — l)th and the (N — k)th 
orbital (see Fig. 11 as an example); these orbitals are 
occupied by spin up electrons only in state jC^y^, x ). 
The non-compact ground states resemble HF states with 
simultaneous charge and spin reconstructions. 9 



B. Separation between bulk and edge electrons and 
exact diagonalization study of large systems 

Through our exact diagonalization study for small sys- 
tems in v < 2 regime, we realize that with smoother 



confinement the compact states with increasing total 
spin become ground states in some tuning parameter re- 
gions; the ground states appear in between compact state 
\Cjif_f.) and \C^~_} k+1 ) have their major charge and spin 
reconstructions in the same orbital segment occupied by 
single spin electrons in compact state j^/v-fc+i)- Assum- 
ing these two points still hold for larger systems, when 
studying the low energy edge excitations and edge insta- 
bilities of v = 2 unpolarized state of large system, we can 
use the exact diagonalization to study a separated edge 
electron system because the occupation number of the 
orbitals inside this edge system (bulk orbitals) are nearly 
2 (completely occupied). The effective Hamiltonian for 
the edge electron system in v < 2 regime is 



mn,aa' C m-\-l,cr C n,cr' c n+l,a' Crn,<r 



m>mo ,n>mo ,1,(7,(7' 

E(ATT cp 4- IJ HF ) 

m>mo ,<7 



C rn,a C m,cr- 



(20) 



The first term is the same as the Coulomb interaction 
term in Eq. (3), except that here we only consider 
the edge electrons in orbitals m > tuq. The orbital 
mo is suitably chosen to study the edge instabilities if 
( n m ) ~ 2 in the exact diagonalization result and the 
choice of mo in our calculation will be stated below. The 
electrons in orbitals m < mo are defined as bulk elec- 
trons in our treatment, mo = (N — N ec i ge )/2, where 
N e dge is the particle number of edge electron system 
we choose and N is the particle number of the whole 
system. AU% mo = U% - U% , where U% is the one 
body confining potential term defined in Sec. II. For 
positive background charge and linear confining poten- 
tial AU^ mo is dependent of particle number N, while 
it is independent of N for parabolic confining potential. 
When edge instability happens, bulk orbitals are still 
completely occupied. So the Coulomb interaction be- 
tween bulk and edge electrons has only Hartree and Fock 
terms left, which is equivalent to a one body Hartree- 
E (2V£, m -^- m ') felt by edge 

m'<mo 

electrons, and V r „, m and V^^ m were defined in Eq. (4). 

Actually the two assumptions at the beginning of this 
subsection for large electron systems are not very secure. 
In disk geometry the sharp peak of the wave function at 
the rath orbital <f) m {r) is at r = \^2mlB, and the wave 
function has width ~ The larger the particle num- 
ber, the more electron orbitals are within the width Ib 
and may be involved in the instabilities in v < 2 regime. 



Fock field U™ 



For the phase transition from unpolarized state 



C 



N/2 
N/2 



to S = 4 compact state 



C 



N/2- 



in our calculation, we 



-4^ 

choose the particle number of edge electron system as 
large as N edge = 18 (i.e. m = N/2 - 9) for N > N edge , 
and give edge electrons 18 orbitals. By this choice, the 
occupation number of the innermost edge orbital (n mo ) 
is larger than 1.99994; the occupation number of the out- 
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FIG. 10: (Color online) Phase diagrams for 10-electron systems under (a) v = 2 positive background charge, (b) parabolic, 
(c) linear confining potential in 1 < v < 2 regime. Step values are A(d/l B ) = 1(T 2 , A6 = 10~ 3 and Aa = 10~ 2 for the three 
confinements respectively. Compact states are labeled by black points, while non-compact states are labeled by red points. The 
total orbital angular momentum number M of each state is labeled above the points. See text for the definition of compact 
state. 
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FIG. 12: (Color online) Low energy spectrum of 30-electron 
system, showing the low energy excitations of v = 2 unpo- 
larized state under parabolic confinement when tuning pa- 
rameter b — 0.14 (mo = 6 in the calculation). The total spin 
quantum numbers S of eigenstates are labeled by different col- 
ors as the annotation. The ground state is v = 2 unpolarized 
state. The dashed boxes enclose the low energy excitations 
in each subspace close to the one of v — 2 unpolarized state, 
and the amount of the states in each box is labeled above it. 
The states pointed by arrows are compact states. 



M 



FIG. 11: (Color online) Occupation numbers (n m ) of (a) com- 
pact state Cf), (b) non-compact ground state (GS) with 
M = 25, S — 1 appearing between Cf) and |C|) in Fig. 
10(b) and (c) compact state |Cf ). Red and blue bars label 
the occupation numbers of spin up and down orbitals respec- 
tively. 



ermost given orbital (n mo +i7) is l ess than 0.00005 for all 
our numerical results up to particle number 80. Because 
bulk electrons do not have tendency to move outward, 
and larger edge orbitals have no tendency to be occupied, 
all the orbitals with occupation numbers not close to or 
2 are included in the edge electron systems we studied. 



Besides, we do not need to worry about the validity of 
the two assumptions for large particle number N up to 
80, since for our choices of mo the exact diagonalization 
results and the assumptions are self-consistent. 

In v < 2 regime, low energy spectra of electron sys- 
tems with particle numbers up to 80 under the three 
confinement are obtained by using exact diagonalization 
on the edge electron systems illustrated above. We will 
use the 30-electron system's low energy spectrum under 
parabolic confinement with b = 0.14 {y = 2 unpolarized 
state is the ground state) to show some common charac- 
teristics shared by these spectra. AM = M — Mq is the 
difference of orbital angular momentum quantum number 
M compared to the one M = Mq for v = 2 unpolarized 
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state (Mo = 210 for N — 30 system). In the subspaces 
with AM > 0, a set of low energy excitations (enclosed 
by dashed boxes in Fig. 12) are separated from higher en- 
ergy states. One of these excitations becomes the ground 
state and thus destabilizes the v — 2 unpolarized state 
with smoother confinement. Similar to the low energy 
excitations of v = 1 ferromagnetic state, we notice that 
the low energy spectrum of v = 2 state is the same as 
the spectrum of bosons excited from the v = 2 state af- 
ter counting the number of low energy excitations in each 
subspace AM > 0. As shown in Fig. 12, from subspace 
with AM = 1, the numbers of low energy excitations are 
1, 2, 3, 5, 7, 11, 15, • • • , the same as the counting of 
bosons. But different from the low energy excitations of 
v = 1 ferromagnetic state, there is no (immediate) simple 
relation between the total spin quantum number S and 
the number of bosons present. Very interestingly, we ob- 
serve that in each subspace M which allows the existence 
of a compact state, the compact state (labeled by arrows 
in Fig. 12) always has the lowest energy. For this rea- 
son, in the following we view the compact states as the 
2nd generation descendants of the v = 2 singlet state, 
and view some low lying non-compact states as their de- 
scendants. This provides us with an essentially complete 
understanding of the low-energy spectrum. 
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FIG. 14: (Color online) Sketch of the orbital configuration of 
compact state 



ICi 1 !) with S z = S. 



FIG. 13: (Color online) Low energy spectrum of 30-electron 
system, showing the low energy excitations on top of com- 
pact ground state | Cig ) under parabolic confinement when 
tuning parameter b = 0.13 (mo = 6 in the calculation). The 
different total spin quantum numbers S of the eigenstates are 
labeled by different colors as the annotation. The ground 
state is compact state |Cif ). The dashed boxes enclose the 



and 



low energy excitations of compact ground state |Cif) 
the amount of the states in each box is labeled below it. 



As it turns out, in each subspace M, part of the low 
energy excitations can be mapped onto bosonic ESW ex- 
citations on top of a neighboring compact state, and just 
like in the v — 1 ferromegnetic state, each boson re- 
duces the compact state's total spin by 1. Taking the low 
energy spectrum of 30-electron system under parabolic 
confinement around the subspace of compact state (C^f ) 
with M = Mo = 219 (Fig. [TBI as an example, in sub- 



spaces with AM = M — Mo > some of the low en- 
ergy excitations are mapped onto the AS = — 1 bosons 
with positive angular momentum q on top of compact 
state jC^f ) (labeled by dashed boxes); in subspaces with 
AM < some of the low energy excitations are mapped 
onto the AS* = — 1 bosons with negative angular momen- 
tum q on top of compact state jC^f ) (labeled by dotted 
boxes). These compact state's bosonic excitations may 
mix with other low energy states in the same subspace, 
which are the excitations of other states in the compact 
state's subspace. The low energy excitations of a com- 
pact state are similar to those of v = 1 ferromagnetic 
state, because both of them can be mapped onto the ex- 
cited bosons which reduce the system's total spin by one 
on top of the original state. This is not surprising since 
a compact state can be viewed as a v — 2 unpolarized 
state surrounded by a v = 1 ferromagnetic state at the 
boundary as shown in Fig. 14. In AM > subspaces, 
excitations of a compact state are determined by the ex- 
citations of bulk v = 2 state and edge ferromagnetic 
state and also their correlation. In AM < subspaces, 
the bulk v — 2 state is inert and excitations of a com- 
pact state are only determined by the edge ferromagnetic 
state. In either case the ferromagnetic state at the edge 
of the compact state plays an essential role on the exci- 
tation. By our exact diagonalization result, we found the 
excited bosons on top of the 5 ^ compact states close 
to v = 2 state have minimum at q = 1 for positive q and 
at q = —1 for negative q. Therefore, the instabilities of 
M a compact state close to v = 2 are just the condensation 
of AS = — 1 , q — 1 bosons (ESWs) at the smoother con- 
finement side and the condensation of AS* = —1, q = —1 
bosons (ESWs) at the stronger confinement side. Two 
neighboring compact states \C^_ k ) and |C^ r ~ 1 fe+1 ) can 
be connected by their excited bosons. As shown in Fig. 
13, the lowest energy state in M = 222 subspace with 
S = has three excited q = 1, AS = — 1 bosons on 
top of compact state |C||); while it can also be viewed 
having four excited q = — 1, AS = — 1 bosons on top of 
compact state |C|g). The lowest energy excited boson of 
compact states close to v — 2 having q = ±1, explains 
our finding in last subsection that the ground states ap- 
pear in between compact state \C^_ k ) and \C^~_} k+1 ) 
have their major charge and spin reconstructions on top 
of these compact states in a segment of orbitals which are 
only occupied by spin up electrons (the orbitals of edge 
ferromagnetic state) in state \C^~_} k+1 ). 

The phase diagrams in v < 2 regime are obtained from 
exact diagonalization for the particle number up to 80 un- 
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FIG. 15: (Color online) Phase diagrams in v < 2 regime for particle numbers from 10 to 80 (even numbers) under (a) v = 2 
positive background charge, (b) parabolic and (c) linear confining potential. Edge electron number N E d ge is chosen to 18 for 
N > 18, and N e d ge = N for N < 18. The orbital number for edge electrons is 18, so electrons can occupy from the moth 
to the (mo + 17)th orbital of the whole system. Step values are A{d/ls) = 10 -3 , A& = 10 -5 and Aa = 10~ 5 for the three 
confinements respectively. Different kinds of states are labeled by the bars with different colors. In the annotation the bosons 
are the ones excited from the compact state in the stronger confinement side of the phase diagram if q = 1 and the compact 
state in the smoother confinement side if q = — 1. Singlet state can be view as the excited bosons on top of the compact states 
in both sides. 



der the three types of confinement potentials in Fig. 15. 
They are qualitatively similar to those of small systems 
in Sec. IV. A. Compact states (gray regions in Fig. 15) 
with different total spins become ground states in some 
tuning paramter regions, and between two neighboring 
compact states, their bosonic excitations (ESWs) with 
q = 1 or — 1 may also destabilize them and appear as 



ground states. Singlet states which can be viewed as the 
condensation of maximum number of bosons from either 
of the nieghboring compact states which always appear 
as ground states (brown regions in Fig. 15). The bosonic 
excitations (ESWs) of compact states with 5^0 may or 
may not destabilize the compact states to become ground 
states. The windows of these compact states' ESWs in 
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phase diagrams will change with different particle num- 
bers, and different confinement types with the step values 
of our calculation. Up to 80 electrons, our exact diago- 
nalization result shows that the initial v — 2 instability 
is toward the S = 1 compact state. As shown in Fig. 15, 
v = 2 edge instability critical point under positive back- 
ground charge confining potential approaches a constant 
d c /ls ~ 0.223, which indicates the existence of ther- 
modynamic limit predicted by our electrostatic model. 
The v = 2 instability critical point b c (a c ) for parabolic 
(linear) confinement decreases (increases) with increas- 
ing particle number, and these are also consistent with 
the prediction of the electrostatic model. 
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FIG. 16: Particle number N vs renormalized critical g factor 
g c — g c fiBB/(e 2 /cIb) in which the edge spin wave (ESW) 
phase window between neighboring compact states closes un- 
der positive background charge confinement. Circular points 
label the close of ESW phase window between S = 1 and 
S = 2 compact state; rectangular points label the close of 
ESW phase window between S — 2 and S = 3 compact state; 
triangular points label the close of ESW phase window be- 
tween S = 3 and S — 4 compact state. 

In the phase diagrams Figs. Qj)l the ESW ground 
states have smaller spin quantum number S compared 
to their two neighboring compact ground states. With 
finite Zeeman coupling, the ESW phase window will 
close and only compact states appear in the phase dia- 
grams. We calculated the renormalized critical g factors 



9c = .9cMs-B/(e 2 /eZs) in which the ESW phase windows 
between neighboring compact states close for different 
particle numbers N under positive background charge 
confinement (Fig. [TBI) . g c approaches a constant ~ 10~ 4 
with increasing particle number. Therefore we predict 
that in thermodynamic limit and under positive back- 
ground charge confinement, with a small Zeeman cou- 
pling (g c ~ 10~ 4 ), the windows of ESW phase close and 
only compact states appear in the phase diagram near 
v = 1, 

V. CONCLUDING REMARKS 

In this paper, we investigate the low energy spin exci- 
tations and edge instabilities triggered by them in integer 
quantum Hall liquids, and the resultant edge reconstruc- 
tions. We conclude that there are likely very rich spin 
structure for the edges of both v = 1 and v — 2 quan- 
tum Hall liquids. We also find that the specific form of 
the instability and the resultant edge structure is very 
sensitive to the details of the confining potential. 

While we have the real electron spin in mind in this 
work, our study and likely many of our results can be gen- 
eralized to systems with pseudospins, including bilayers 
quantum Hall systems 1 ^ and graphene*2P In some sense 
these systems can be even more interesting. For exam- 
ples there is no Zeeman coupling to the pseudospins and 
thus no associated energy penalty for pseudospin tex- 
tures; the Dirac Landau level wave function allows for 
spin and pseudospin textures in higher Landau levels^ 
and the lack of SU(2) symmetry allows more exotic forms 
of pseudospin texture in bilayer systems^ 

Another natural direction to pursue is edge spin exci- 
tations and textures in fractional quantum Hall liquid. 
In fact a preliminary attempt has already been made in 
that direction^ 
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